Quantum-tunneling dynamics of a spin-polarized Fermi gas in a double-well potential 
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We study the exact dynamics of a one-dimensional spin-polarized gas of fermions in a double- well 
potential at zero and finite temperature. Despite the system is made of non-interacting fermions, 
its dynamics can be quite complex, showing strongly aperiodic spatio-temporal patterns during 
the tunneling. The extension of these results to the case of mixtures of spin-polarized fermions in 
interaction with self-trapped Bose-Einstein condensates (BECs) at zero temperature is considered 
as well. In this case we show that the fermionic dynamics remains qualitatively similar to the one 
observed in absence of BEC but with the Rabi frequencies of fermionic excited states explicitly 
depending on the number of bosons and on the boson-fermion interaction strength. From this, the 
possibility to control quantum fermionic dynamics by means of Feshbach resonances is suggested. 

PACS numbers: 03.70.Lm; 03.75.+k; 03.05.Ss 



I. INTRODUCTION 

The control and manipulation of confined ultracold 
atomic gases make possible the study of the dynamics 
of many-body quantum states either for bosons [l|-|5| and 
for fermions 0,0], as well as for their mixtures. Ultracold 
atoms give the opportunity to investigate the most im- 
portant quantum effects with a very high level of control 
and precision [8l-[ll|. Among these effects, the quantum 
tunneling is a topic of wide interest. In this sense, one of 
the most important examples is provided by Bose gases 
confined in double-well shaped potentials. In these sys- 
tems the tunneling of bosons through the central barrier, 
which separates the two potential wells, is responsible 
for the atomic counterpart of the Josephson effect (see, 
for instance, Refs. jll| - |l4| ) observed in superconductive 
junctions |15| . The macroscopic quantum tunneling was 
studied also with two weakly-linked superfluids made of 
interacting fermionic atoms for which it is possible to 
obtain atomic Josephson junction equations |16| . 

The optical lattices and double- well traps are paradig- 
matic external potentials to understand the tunneling in 
atomic systems characterized by a small number of parti- 
cles [l|, \4. \it\ . The dynamics of small systems of spinlcss 
bosons llSl. of finite Fermi-Hubbard and Bose-Hubbard 
systems [l9[, and of spin 1/2 fermions - especially within 
the quantum information community [20l . [2]| - were 
widely studied. 

In this work we consider a one-dimensional (ID) dilute 
and ultracold gas of spin-polarized fermionic atoms in a 
double-well potential in the absence and in the presence 
of a Bose-Einstein condensate (BEC) which is intrinsi- 
cally localized in one of the two potential wells. In both 
cases we have that in correspondence to certain values of 
the parameters of the double- well potential, its eigenspec- 
trum is organized in doublets spaced between themselves 
and made of quasi-degenerate states. Each doublet is 
characterized by its own Rabi frequency, proportional to 
the gap between the two corresponding states. It turns 



out that even in absence of BEC and of interactions be- 
tween the fermions, this kind of system exhibits an inter- 
esting behavior. 

In the first part of the paper we focus on a spin- 
polarized quantum gas in absence of BEC. We suppose 
that at the initial time the system is prepared in such 
a way that all the fermions are localized at a given side 
of the barrier and that a number of fermionic states can 
be excited by means of external resonant fields (e.g. res- 
onant with the Rabi frequencies of the doublets). The 
dynamics sustained by the above Hamiltonian allows to 
study the transfer of matter which takes place through 
the barrier by tunneling. This effect was studied in 
Ref. [l8| for few bosons in a ID double-well at zero- 
temperature. Indeed, the non-interacting pure fermionic 
system can be viewed as a realization in the infinitely 
strong repulsive interaction limit of ID hard-core bosonic 
systems j22| . where the Pauli exclusion principle mimics 
the hard-core interaction. By following Ref. [l^l, we fix 
the number of fermions and analyze the spatial variations 
of the density profile of the fermionic cloud at different 
times, and the total fermionic fractional imbalance Zp 
between the two wells as a function of time t. We carry 
out this analysis both with one fcrmion and with more 
fermions by pointing out the striking differences between 
the two cases. In the former case the temporal behavior 
of ZF{t) is fully periodic and its period (Rabi period) is 
given by the inverse Rabi frequency of the lowest dou- 
blet. In the latter case the periodicity occurs on much 
longer times, equal to the minimum common multiple of 
the Rabi periods of the populated doublets. Zp{t) may 
even be completely aperiodic if some of the Rabi periods 
of the populated doublets arc incommensurate. In any 
case, the shape of the density profile generally exhibits 
spatio temporal patterns which do not replicate the ini- 
tial situation for rather long times. In Refs. [l^ fisj - 
the atomic systems therein considered are dealt with at 
zero-temperature. Here we perform our analysis also at 
nonzero temperature. By keeping fixed the number of 
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fermions, we study the influence of tlie temperature on 
the total fractional imbalance: there is a spatial broad- 
enings of the density profile that become more evident 
by increasing the temperature. 

In the second part of the paper we discuss the case of 
a spin-polarized quantum Fermi gas interacting with a 
quasi-stationary BEC at zero temperature. In particu- 
lar, we show that for a self-trapped BEC (e.g. a BEC 
intrinsically localized by the nonlinearity in one of the 
two wells of the potential) and for small boson-fermion 
interactions, the fermionic imbalance dynamics remains 
qualitatively similar to the one of the pure fermionic 
case with the only difference that the spacings of the 
fermionic levels (and therefore the Rabi frequencies) be- 
come dependent on both the number of bosons and the 
boson-fermion interaction. This suggests the possibility 
to control the fermionic quantum dynamics by changing, 
for example, the inter-species scattering lengths by using 
Feshbach resonances, a fact which could be of interest for 
applications to quantum computing. 



II. THE SYSTEM 



where (j)aj{x) are the complete set of ort ho normal eigen- 
functions and j the corresponding eigenenergies. Here 
a = 1, 2, 3, ... gives the ordering number of the doublets 
of quasi-degenerate states, and j = S,A gives the sym- 
metry of the eigenfunctions {S means symmetric and A 
means anti-symmetric). It follows that the eigenvalues 
are ordered in the following way: ei_s < ei^A < £2.s < 
^2.A < £3,5 < £3, A < ■■■ ■ Without loss of generality, we 
consider real eigenfunctions (t>a.iix). 

Due to the symmetry of the problem, it is useful to 
introduce a complete set of orthonormal functions 



1 



(4) 



and 



£.a,Lix) = --j={(j)a,s{x) ~ (l)a,A{x)) . (5) 



a.A, such that 



If we fix the phase of (j)a,s ^nd 
(/)^ j(-t-oo) < for both j ^ S and j = A then S^^x 
and S,a,R are mainly localized in the left and right well 
respectively, at least for the low lying states. 

The field operator ip{x) can be written as 



We consider a confined dilute and ultracold spin- 
polarized gas of Np fermionic atoms of mass m in a 
double- well potential. For dilute and ultracold atoms the 
only active channel of the inter-atomic potential is the s- 
wave scattering. For spin-polarized fermions the s-wave 
channel is inhibited by statistics and consequently the 
system is practically an ideal Fermi gas. The trapping 
potential of the system is given by 



Q=l i=LM 



(6) 



in terms of the single-particle fermi operators Ca^i (c]^ j) 
satisfying the usual anticommutation rules, which de- 
stroy (create) one fermion in the state ^a.L{x) or £,a.B,i^)- 
By inserting Eq. © into ([5]), the Hamiltonian takes 
the form 



V{r) = Uix) + ljnu;l{y^ + z^) , (1) 

We suppose that the system is one-dimensional (ID), 
due to a strong radial transverse harmonic confinement 
[23I [23 | , with a symmetric double- well trapping potential 
in the longitudinal axial direction, that we will denote by 
U (x) . We are thus assuming that the transverse energy 
hLj± of the confinement is much larger than the charac- 
teristic energy of fermions in the axial direction [2^, . 

On the basis of the dimensional reduction, the ID 
Hamiltonian of the dilute spin-polarized Fermi gas con- 
fined by the potential U{x) is 



H 



4;\x) 



2m dx^ 



U{x) 



i){x) dx , (2) 



where '<p{x) and V'^(a;) are the usual field operators an- 
nihilating or creating a fermion at position a;, therefore 
obeying anticommutation rules. The single-particle sta- 
tionary Schrodinger equation associated to the Hamilto- 
nian ^ can be written as 



(l^ajix) = ta,]4>a,]{,x) , (3) 



Q — 1 i—L,B. a— 1 

(7) 

where hF,a,i = c^a fermionic number operator 

of the a-th state in the i well {i = L,R: L means left 
and R means right). Notice that the diagonal part of the 
Hamiltonian ([7|), given by 



(8) 



a=l i=L,R 



describes a double- well Fermi system with an infinite bar- 
rier, i.e. with zero hopping terms. The energy 



given by 



£,a,i {x) 



2m dx^ 



U{x) 



(9) 



(10) 



is the same for the left and right states described by 
^a,Lix) and ^a.ni^) The energy 



(11) 
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is the energy of hopping from the L to the R well within 
the same doublet and gives directly half of the Rabi fre- 
quency of the a-th doublet as 



fin 



J a 



(12) 



i.e. the Rabi angular frequency 51^°'" is given by fi^"''* 
2fia = 2Ja/h [2|. 



with < /q^l < 1 and Q < fa,R < ^ parameters fixed by 
the choice of the statistical density operator. We remark 
that fa,L and /q, ^ can be any distribution. In particular 
they could be the initial distributions corresponding to 
thermal equilibrium in the two fully separate wells with 
different number of particles in each well. In this case 
the statistical density operator is 



Po = e 



{Ho-tJ.F,LNF,L-p,F,RNF_R) 



(22) 



III. TUNNELING DYNAMICS OF FREE 
FERMIONS IN DOUBLE WELL POTENTIAL 

To study the dynamics sustained by the Hamiltonian 
([7|), we analyze the problem within the Heisenberg pic- 
ture. The Heisenberg equation of motion for the operator 
ip{x,t) is 



^n—'^p[x,t) 



2m dx^ 



U{x) 



and the time-dependent density operator is 
fip{x,t) = ?/'^(a;, <)-0(a;, i) . 



■0(x,i) , 
(13) 

(14) 



The Heisenberg equations of motion for the operators 
Tip a L and fiFa.R rcad 



ifl — flFa.L = -Ja{ha.L 

at 



ih — nF,a.R 

at 

at 



(15) 
(16) 



^a,R 



Ca,L- The equa- 



where h^^L = c^a.L^a,R^ ^^^^ ha^R 
tions of motion for the operators hax and ha^R are 



ifl — ha^L — —Ja{nF.a,L—nF,a,R) 

at 



iTT'-Trha.R 

at 



Ja{flF,a,L — flF,a,R) 



(17) 
(18) 



We want to analyze the spatio-temporal evolution of 
our system both at zero and finite temperature T. To 
this end, we evaluate the thermal average of both sides 
of the Heisenberg equations of motion obtained so far. 
The thermal average of an operator A{t) is given by 



{A{t)) = 



Tr 


i(t)po 


Tr 







= A{t) , 



(19) 



where po is the chosen statistical density operator. It is 
clear that a thermal average evaluated by using the den- 
sity operator p = e-'^(^-^'^) (with P = l/iksT)) does 
not give rise to dynamics (time-dependent observables) 
[26| . We want a statistical operator po which produces 
the initial conditions 



"-F,q,l(0) = (flF.a.LiO)) = fax , 



(20) 
(21) 



where Hq is given by Eq. dS]), TVf,^ = ^^F,Q,^ and 

fiF,i the number operator and the chemical potential of 
fermions in the i-th well, respectively. In this way we 
find 



fa,L = 
fa,R — 



1 



g/9 (ec-MF.t) ^ 1 
1 

(Sc-fJ-F.R) -(_ 1 



(23) 
(24) 



The chemical potentials pp.L and pf.r are fixed by the 
number of particles in the left and right wells at time 
zero: 



1 



a=l 
oo 



Nr,l{0) = E 



gPi^a-tJ.F.L) ^ I 
1 

Ql3(ea-fJ,F,R) _|_ I 



(25) 
(26) 



It is easy to show that the time-dependent density pro- 
file nF{x,t) of the Fermi system. 



nF[x,t) = {t/j\x,t)ip{x,t)) 
can be written as 



uf 



(27) 



(28) 



a=l i=L,R 



where nF.a,i{t) ~ {nF,a.i(t)). In addition, from Eqs. 
([T5 |) - p8)) it is straightforward to get the following cou- 
pled ODEs 

nF,a,L + 2fi^ nF,a,L = nF,a.,R , (29) 

nF,aM. + 2f2^ nF,a,R = 2r2^ nF,aX (30) 

with fla given by Eq. (IT^ . 

Notice that the total number nF,a = nF,a,L{t) + 
nF.a,R{t) of fermions in the a-th doublet and the 
fcrmionic hopping number ha{t) — ha. hit) + ha,R{t) are 
both constants of motion. 

It is not difficult to show that the ODEs (gH) and ^ 
with the initial conditions (PO)) and (PT|) have the follow- 
ing solutions 

nF.a,L{t) = fa,L COS^ (^Q*) + faM siu^ (flat) , (31) 
nF.a,R{t) = fa.R COS^ {^at) + fa,L siu^ {^at) ■ (32) 
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FIG. 1: The double well potential ((35]) for A = 5■10"^ 5 = 1, 
and C = 5. For these values of the double-well potential 
parameters, there are 30 energetic levels - corresponding to 15 
doublets - with energy smaller than the height of the barrier. 
Energies are measured in units of fiuj±, lengths in units of 



muj±. 



The population imbalance Zq.(<) within the a-th double 
is 

ZF.ait) = nF,a.L{t)~riF.a.R{i) = ifa.L^fa.R) COS (217Qi) , 

(33) 

and the total ferniionic imbalance ZF{t) is given by 

^^^^^ = ^ E -^"W = 1^ Y.^fo.,L~fo.,R) COs(21)„t) . 



1=1 



(34) 

We model the double- vifell trapping potential U{x) (see 
Fig. 1) in the form 



U{x) = Ax^ + B{e 



1) 



(35) 



We study the time-dependent density profile (P7)) and 
the total fermionic imbalance p4p both at zero and finite 
temperature for different number of fermions. The results 
of this study are shown in Fig. 2 (Np — 1), Fig. 3 
{Nf = 2), Fig. 4 {Nf = 6), and Fig. 5 \Nf = 12). 

In our calculations we have supposed that the fermions 
of the system are atoms, and set the harmonic trans- 
verse frequency w_l to 160 kHz, so that — 0.1/x m j23 |. 
The absolute temperatures used in Figs. 2-5 are calcu- 
lated according to the above assumptions. In each of 
these figures we show, for fixed values of T, the spatial 
evolution of nF{x,t) at different times and the temporal 
evolution of Zf- When Nf = 1 and T = - the continu- 
ous line of Fig. 2 - only the lowest doublet is involved in 
the dynamics. At zero temperature, the temporal behav- 
ior oi Zf is harmonic and is characterized by only one 
frequency: 2fii. At finite temperatures - the dashed and 
the dot-dashed lines of Fig. 2 - the fermions are allowed 
to partially occupy doublets which are empty at T = 
(see the captions of Fig. 2), and the harmonicity of Zf is 
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FIG. 2: (Color online). Density profile np vs. space x and 
total fermionic imbalance Zp vs. time t with Np = 1. Plots 
of np(x) are for different values of t. We define t\ as 27r/f2i. 
In this figure: t = and t = 0.25 1\ (the two upper panels 
of np from left to right); t — 0.88 ti and t = 1.5 (the two 
lower panels of wf from left to right). For both the nF{x) and 
Zp{t) plots, the continuous line represents data for Kp T = 0, 
the dashed line for Kp T — 0.038 (~ (e2 — ei)), and the dot- 
dashed line for KpT = 0.087 (~ (ea - ei)). The plots at zero 
temperature are obtained with the fermion initially localized 
in the ground state of the right well; the plots at finite temper- 
ature are obtained with Np^l{0) = and Np^r{0) = 1 with 
level occupation defined by Eq. (p6|) . Times are measured in 
units of {u)±)~^ , energies in units of hujj_, and lengths in units 
of a± = y^h/m(jj±. 



deformed by this thermal effect, as we can sec from the 
Zf panel of Fig. 2. 

When the number of fermions is greater than one (Figs. 
3-5), doublets other than the lowest one are involved in 
the dynamics. At T = the temporal evolution of the 
total fractional imbalance - the continuous line of the Zp 
plots - is no more characterized by a single frequency as in 
the Np = 1 case; it is, in fact, characterized by a mixing 
of frequencies, each of them proportional to the gap of 
the corresponding filled doublet. As in the Np ~ 1 case, 
at finite temperatures the fermions can partially occupy 
doublets higher in energy than those filled at T = (see 
the captions of Figs. 3-5). The effect is to induce a defor- 
mation in the oscillations of the total fermionic fractional 
imbalance - with respect to T = case - as shown in the 
plots of Zp - the dashed and the dot-dashed lines - of 
Figs. 3-5. The thermal effects influence the fermionic 
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FIG. 3: (Color online). Density profile np vs. space x and 
total fermionic imbalance Zp vs. time t with Np = 2. Plots 
of npix) are for different values of t. We define ti as 27r/Sli. 
In this figure: t = and t = 0.74 ti (the two upper panels 
of np from left to right); t — 0.82 and t = 1.5 ii (the two 
lower panels of np from left to right). For both the np{x) and 
Zp{t) plots, the continuous line represents data for Kb T = 0, 
the dashed line for Kb T = 0.047 (~ {is - £2)), and the dot- 
dashed line for Kb T = 0.164 (~ (es — £2)). The plots at zero 
temperature are obtained with the 2 fermions initially in the 
2 lower states of the right well; the plots at finite temperature 
are obtained with Np^l{0) = and Na^LiO) = 2 with level 
occupation defined by Eq. (|26|l . Units are the same as in Fig. 
2. 
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FIG. 4: (Color online). Density profile np vs. space x and 
total fermionic imbalance Zp vs. time t with Np = 6. Plots 
of np{x) are for different values of t. We define ti as 2n/Qi. 
In this figure: t = and t = 0.76 ti (the two upper panels 
of np from left to right); t — 0.83 and t = 1.52 ti (the two 
lower panels of np from left to right). For both the np{x) and 
Zp{t) plots, the continuous line represents data for Kb T — 0, 
the dashed line for Kb T — 0.068 (~ (ey — es)), and the dot- 
dashed line for Kb T = 0.219 (~ (eg — eg)). The plots at zero 
temperature are obtained with the 6 fermions initially in the 
6 lower states of the right well; the plots at finite temperature 
are obtained with A^f,l(0) = and Npi^l{0) ~ 6 with level 
occupation defined by Eq. (|26|) . Units are the same as in Fig. 
2. 



density profile as well. In particular, when the temper- 
ature is finite, npix) experiences a broadening with re- 
spect to the case T = 0. From the plots oi np shown in 
Figs. 2-5, we see that the higher is the temperature the 
wider is this broadening. 

Finally, it is worth observing that during the dynamics 
of our model the system is not in thermal equilibrium and 
it will never reach it since the particles do not interact 
between themselves. Each doublet preserves its energy 
and its number of particles. 



IV. TUNNELING DYNAMICS OF 
SPIN-POLARIZED FERMIONS INTERACTING 
WITH INTRINSICALLY LOCALIZED BEG 

In this section we consider the case of a spin-polarized 
quantum Fermi gas in interaction with a BEC which is 
intrinsically localized in one of the two wells of the poten- 



tial ([55]) at zero-temperature. In particular, we are inter- 
ested in the changes in the Rabi frequencies and in the 
fermionic imbalance dynamics induced by the presence of 
the BEC. To this regard we restrict, for simplicity, to the 
case of a small number of excited fermionic states present 
in the fermionic density and to weak boson-fermion inter- 
actions. In this situation the condensate will remain self- 
trapped (practically stationary) in the course of time and 
the fermionic imbalance dynamics will be qualitatively 
similar to the one discussed in the previous section. The 
fermionic spectrum and the Rabi frequencies, however, 
will depend on the boson-fermion interaction due to the 
presence of a bosonic effective potential in the fermionic 
Schrodingcr equation of motion (see below). 

To describe the spin-polarized fermionic gas in inter- 
action with the BEC we adopt a mean field description 
for the condensate but treat the fermions still quantum 
mechanically. In this case one can show [2^, [2^ that the 
dynamics of the mixture is described by the following set 
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FIG. 5: (Color online). Density profile np vs. space x and 
total fermionic imbalance Zp vs. time t with Np = 12. Plots 
of npix) are for different values of t. We define ti as 27r/Sli. 
In this figure: t = and t = 0.77 ii (the two upper panels 
of np from left to right); t = 0.83 ti and t = 1.52 ti (the two 
lower panels of np from left to right). For both the np{x) and 
Zp{t) plots, the continuous line represents data for Kb T — 0, 
the dashed line for KbT = 0.142, ((eis - £12) < KbT < 
(ei4 — £12)), and the dot-dashed line for Kb T = 0.333, ((eis — 
£12) < Kb T < (ei6 — £12)). The plots at zero temperature are 
obtained with the 12 fermions initially in the 12 lower states 
of the right well; the plots at finite temperature are obtained 
with Np^l{0) = and Nu^l(0) = 12 with level occupation 
defined by Eq. (|26|) . Units are the same as in Fig. 2. 



of coupled equations 



ih 



9* _ 
'dt ~ 



ih 



dt 



2m dx^ 



2m dx^ 



U{x) + geNeM^ + gBFnp 



-^U{x) + NBgBF\^ 



Xj 



(36) 
(37) 



where np{x,t) = X^j^Ti IXj(2^i^)P denotes the fermionic 
density with Xj{x,t) (j = l,...,Np) the set of orthonor- 
mal wave functions which satisfy Eq. (P7|) . ^(x, t) is the 
bosonic wavefunction normalized to one and such that 
nB{x,t) = NB\'^{x,t)\'^ is the bosonic density with Nb 
the total number of bosons. The ID interaction strengths 
are gB = 2h'^aB/{'ma\) and gBP = 2fr'aBF/{'ma\)^ 
with aB and a^F the boson-boson and boson- fermion 
s-wave scattering lengths and a± = ^Jh/{mLu^) the char- 
acteristic length of the strong transverse harmonic con- 
finement of frequency a;^ [zsl . [2^ . A similar study for 



tion with frozen fermions was done in |27| . For sim- 
plicity, we also assumed equal masses and the same 
trapping potential for both bosons and fermions. The 
bosonic cloud is self-trapped in one of the two wells 
when the nonlinear strength .gsiVs exceeds a certain fi- 
nite critical value [T2j . This condition can be written 
as \aB\NB/ai^ > {ei^a - ei,s)/{huj±) [H], and it is fully 
satisfied in our numerical experiments. 

Note that the BEC wavefunction evolves according to 
a nonlinear Gross-Pitaevskii equation (GPE) in which 
the fermionic density enters as a potential, while the 
fermionic dynamics is still linear with the condensate 
density playing the role of an additional external poten- 
tial. From this it is clear that the fermionic eigenstates 
and eigenvalues can be well approximated by the linear 
Schrodinger equation 



' dt 



2m dx^ 



Xj 



(38) 



with the effective potential Ueff{x) = U{x) + gBF^iBix), 
with nB{x) denoting the stationary bosonic density. We 
have numerically verified that indeed the bosonic cloud 
is practically stationary. From Eq. the possibility 

to manipulate the fermionic dynamics by changing Nb 
or the inter-species scattering length becomes evident. 

In order to check these predictions, we have numer- 
ically determined the stationary states of the mixture 
by solving in a self-consistent manner the time indepen- 
dent equations corresponding to Eqs. ([36]) and ([37)) f29j 
for the case of attractive interactions. In Fig. [S^) we 
show the stationary bosonic wavefunction localized in 
the left well of the potential and the first 30 fermionic 
stationary levels obtained in presence of BEC with the 
self-consistent method. Notice that the lowest levels de- 
viate from the corresponding ones obtained in Fig. [T]and 
the quasi degeneracy of the lowest levels is removed due 
to the bosonic effective potential (in the present case the 
Rabi frequencies of lowest levels are increased due to the 
level splitting). In panel b) of Fig. [5] we show the first ten 
fermionic stationary wavefunctions corresponding to the 
lowest ten energy levels of panel a), while in panel c) we 
depict the first three excited (e.g. non stationary) states, 
C,i{x), constructed from the lowest stationary wavefunc- 
tions as: Q{x) = {x2i-i{x) -\- X2t{x)) / V2. Note that the 
fermionic lowest energy eigenstate is asymmetric, having 
the same localization of the condensate due to the attrac- 
tive Bose-Fermi interaction, while the next level has the 
opposite localization. This implies that the lowest ex- 
cited states are not localized in the same potential well, 
as for the pure fermionic case considered before, but are 
extended between the two wells. 

Excited fermionic densities corresponding to ten 
fermions and including the lowest excited states are 
shown in panel d) of Fig. [6] together with the stationary 
density. These excited densities were used, together with 
the BEC wavefunction shown in the top part of panel 
6a), as initial conditions to calculate the total fermionic 
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FIG. 6: Panel a): condensate localized wavefunction 
(top curve) corresponding to Nb = 470 bosons in a dou- 
ble well potential and in the presence of 30 spin polarized 
fermions with attractive boson-boson and boson-fermion in- 
teractions as = —0.001, fflsF = —0.001. Horizontal lines de- 
note the first 30 fermionic energy levels in the presence of the 
condensate while the bottom curve represents the trapping 
potential with parameters fixed as in Fig. [1] Panel b): lowest 
ten fermionic eigenfunctions Xi{^) (from bottom to top) in 
the presence (continuous curves) and in the absence (dotted 
lines) of the BEG wavefunction depicted in panel a) (an offset 
of 0.5 between curves has been added to avoid overlapping). 
Panel c): first three excited states in the presence of BEG 
constructed with the lowest eigenstates of panel b) (see text). 
An offset of 0.7 between curves was added to avoid overlap- 
ping. Panel d): fermionic densities for 10 fermions including 
one (dotted line), two (dash dotted line), and three (dashed 
line) excited states depicted in panel c. The continuous line 
denote the stationary fermionic density in absence of excited 
states. Energies are measured in units of hu}± , lengths in units 
of a± — y^h/mij±. The bosonic and fermionic wavefunctions 
are both normalized to one. 



imbalance density, Zp{t), from direct nmnerical integra- 
tions of Eqs. and ([57)) . The results are depicted 
in Fig. [3 panels a)-c). From panel a) of this figure we 
see that when only one excited state is present in the 
density, the imbalance dynamics is periodic with a single 
frequency in the spectrum. The Rabi linear frequency 
i^a ~ O,^"-^^ / {2tt) = ila/n calculated from the dynamics 
of the fermionic imbalance density is found in very good 
agreement (see panel d)) with the value obtained from 
the energy spectrum and calculated from the nonlinear 
eigenvalue problem associated with Eqs. (|36|) and ([37)) 
using the self-consistent approach [2^. The same good 
agreement is found in the case of two and three excited 
states present in the fermionic density panels (compare 
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FIG. 7: Panels a)-c): dynamics of the fermionic density im- 
balance Zpii) (bottom) and corresponding Fourier spectrum 
(top part of panels) as obtained from numerical integrations 
of Eqs. (|36|) and (|37|l using as initial conditions the station- 
ary localized bosonic wavefunction and the non stationary 
fermionic densities with one (a), two (b) and three (c) excited 
states depicted in Fig. [6] Panel d): dependence of the Rabi 
linear frequency Va on the boson-fermion interaction for the 
first five excited levels indicated with dots, squares, plus, tri- 
angle and crosses, respectively. The vertical line indicates the 
value of aBF at which the other panels have been evaluated. 
Parameters are fixed as in Fig. |6] Time is measured in units 
of ctJ_L~^ and asF is measured in units of a^. 



panels b),c), with corresponding values of panel d)). 

From panels a)-c) of Fig. |7]it is clear that the fermionic 
imbalance dynamics is given by a superposition of har- 
monics in number equal to the number of excited states 
and with periods fixed by the level spacing Rabi frequen- 
cies. Also note from panel c) of Fig. [7] the presence of 
beatings of period l/(i^2 — ^3) generated by the the two 
close frequencies 1/2 and , as expected for linear systems 
(notice also the presence of small peaks close to the Rabi 
brequcncy ui probably due to the quasi-stationarity of 
the bosons). 

The dependence of the Rabi frequency on the Bose- 
Fermi interaction is depicted in Fig|7)i) for the first five 
excited states. Note that the frequency of the first excited 
state has a large variation with asF due to the fact that 
the corresponding eigenstates are the ones most effected 
by the presence of the condensate (see Fig. [6)3). Similar 
dependences are also obtained by changing the number 
of atoms in the condensate and keeping fixed the fixed 
boson-fermion interaction. 

From this we conclude that in the BEC self-trapped 
regime the dynamics for the fermionic density imbalance 
remains qualitatively similar to the one observed in the 
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pure fermionic case but with the Rabi frequencies explic- 
itly dependent on the boson-fermion interaction. Since 
this interaction can be easily changed by external mag- 
netic fields via Feshbach resonances, we have found an 
effective way to control the the dynamics of the spin po- 
larized fermionic gas which can be implemented in a real 
experiment. 

V. CONCLUSIONS 

We have investigated a confined dilute and ultracold 
spin-polarized gas of fermionic atoms in a double-well 
potential. We have analyzed the quantum tunneling 
through the central barrier by studying the density pro- 
file and the total fermionic fractional imbalance. We have 
pointed out that despite the fermions do not interact be- 
tween themselves, the dynamics is quite complex, and it 
exhibits very interesting aperiodicity features. We have 
performed our analysis both at zero and finite tempera- 
ture by discussing the role of the temperature in produc- 



ing deformations of the density profile with respect to the 
zero-tempcrature case. We also have discussed the pos- 
sibility to include in the system a bosonic component, 
which, under given hypothesis, does not produce impor- 
tant changes in the dynamics of the system at least within 
very low temperature regimes. In particular, we have 
shown that the presence of a self trapped Bose-Einstein 
condensate weakly interacting with the spin-polarized 
fermi gas allows to achieve a quite effective control of the 
Rabi frequencies of excited fermionic states by chang- 
ing the boson-fermion interaction with external mag- 
netic fields via Feshbach resonances. This Bose-Einstein 
condensate-induced control of the fermionic quantum dy- 
namics and can be implemented in a real experiment and 
can be of interest for applications to quantum comput- 
ing. 
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